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1 Short-time dynamic scaling 

For a critical dynamic system in the long-time regime, it is well known that there 
exists a universal dynamic scaling form. This is more or less due to that the 
correlation time is divergent and the spatial correlation length is also very large 
in the long-time regime. 

Is there any universal behavior in the macroscopic short-time regime of the 
critical dynamic evolution? The traditional answer is no. However, this has 
been changed in recent years. In the macroscopic short-time regime, the spatial 
correlation length is also small in the macroscopic sense. However, the large 
correlation time induces a memory effect. The memory effect is represented by a 
scaling form. 

Let us consider that a magnetic system initially at high temperature and with 
a small magnetization is suddenly quenched to the critical temperature, then re- 
leased to dynamic evolution of model A. A universal dynamic scaling form, which 
sets in right after a time scale tmic which is enough large in microscopic sense but 
still very small in macroscopic sense, has been derived with an e-expansion |l[]. 
Important is that extra critical exponents should be introduced to describe the 
dependence of the scaling behavior on the initial conditions. For example, for the 
k-th moment of the magnetization, the finite size scaling form is written as |I[] 

M('=)(t, r, L, mo) = h-^f^'" M^^\h-H, h^'^r, b'^L, b'^'mo). (1) 

Here /3, v are the well known static critical exponents and z is the dynamic 
exponent, while the new independent exponent Xq is the scaling dimension of the 
initial magnetization m^. 

In this and the next sections we discuss the dynamics generated by Monte 
Carlo algorithms, which belongs to model A ^j. 

For a large lattice (L = cxd) and r = , from the scaling form (|l]) one derives 
for small enough t and mo fl], |^ 

M(t)~mot^ e={xo-(3/u)/z. (2) 

For almost all statistical systems studied up to now, the exponent 9 is positive, 
i.e. the magnetization undergoes surprisingly a critical initial increase [|l|, |^. The 
time scale of this increase is to ~ m,o^ °. In Fig. |l] (a), the magnetization for 
the 3D Ising model with the heat-bath algorithm has been displayed in log-log 
scale. Power law behavior is observed at beginning of the time evolution. The 
microscopic time scale tmjc is extremely small in this case. Extrapolating the 
results to mo = one obtains 9 = 0.108(2) §. 

Other two interesting observables in short-time dynamics are the auto-correlation 
and the second moment of the magnetization. For r = and mo = 0, it is well 
known d, I 

M(2)(t)~t-2^ C2 = id-2l3/iy)/z. (3) 

1 



M(t) 



0.1 



m„=o.06 






^^^^.-^--^ 0.04 


— 


-— " 


0.02 


_- 


-^ 



10"' 



M'"(t) 



10 



10"' 




10 



100 



10 



100 



Figure 1: (a) The magnetization of the 3D Ising model with L = 128. (b) The 
second moment. 



d 



Ca 



C2 



2(3/v 



0.108(2) 1.36(1) 2.04(2) 0.970(11) 1.02(3) 
Table 1: The exponents of the 3D Ising model. 



In Fig. |l] (b), the second moment for the 3D Ising model with the heat-bath 
algorithm is plotted in log-log scale [Q. Almost perfect power law is seen. Careful 
analysis reveals the auto-correlation (r = 0) 



A{t) ~ r 



z 



(4) 



Interesting here is that even though vtiq = 0, the exponent 6 still enters the 
auto-correlation. The behavior in Eq. (^) has been confirmed in many systems 

ii- 

In Table |l|, we have summarized the measured critical exponents of the 3D 

Ising model. Taking the exponent 6 as an input, from Ca we estimate the dynamic 

exponent z. With z in hand, from C2 we obtain the static exponent 2(3 /v. The 

values of z and f3/v agree nicely with those measured in equilibrium, z = 2.04(3) 

[0 and 2/3/i/ = 1.036(14) 0. This fact, on the one hand, strongly supports the 

short-time dynamic scaling, and on the other hand, provides new ways for the 

numerical measurements of the critical exponents. 



2 Dynamic measurements of critical exponents 

For critical dynamics, correlation times are very large. Therefore, for standard 
Monte Carlo simulations in equilibrium it is very difficult to generate independent 
spin configurations. This is the so-called critical slowing down. One of the most 
succeesful methods to overcome critical slowing down is the non-local cluster 
method. However, it does not apply to systems with quenched randomness or/and 
frustration. 

In short-time dynamics, we do not have the problem of generating independent 
spin configurations. Therefore, one does not suffer from critical slowing down. 
The 2D fully frustrated XY (FFXY) model is a typical example where critical 
slowing down is severe. 

In general, for determination of the dynamic exponent z and static exponents a 
dynamic process starting from a completely ordered state is more favorable, since 
fluctuation is much less. The completely ordered state in the 2D FFXY model 
is the ground state where spins on four sublattices orient in different directions. 
Here we only concern the chiral degrees of freedom and the phase transition is 
of the second order. Assuming the lattice is sufficiently large, around the critical 
point the dynamic scaling form of the magnetization is written as 

M{t, t) = t-'^I^^F{t^l''^T) . (5) 

If r = 0, the magnetization decays by a power law M{t) ~ t^f^/'^^ . If r 7^ 0, the 
power law behavior is modified by the scaling function F{t}^'^^T). From this fact, 
one determines the critical point and the critical exponent (3/vz. 

In Fig. 1^ (a), the magnetization at T = 0.452, 0.454 and 0.456 is plotted 
in double-log scale. The Metropolis algorithm has been used in simulations. 
Extrapolating M{t) to other temperatures, from the time interval [200,2000] we 
measure the critical temperature to be Tc = 0.4545(2) and P/i'z = 0.0602(2) 0. 
To determine the exponent l/{i>z), we differentiate In M{t,T) and obtain 

dr In M(t,r)|.=o~t'/'^'. (6) 

This power law behavior is shown in Fig. (b) 0. 

In order to estimate the dynamic exponent z independently, we introduce a 
Binder cumulant U = M'^^-'/M^ — 1, and finite size scaling analysis shows 

U{t,L)r~.t^/'. (7) 

Now we complete the measurements of the critical exponents. The results are 
given in Table H, in comparison to those obtained in equilibrium simulations. 
Our results obviously support that the exponent u of the 2D FFXY model is 
different from that of the 2D Ising model. But other exponents of both models 
look similar. 
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0.191(3) 



Table 2: Dynamic measurements of the exponents of the 2D FFXY model 
compared with those in recent references. For the Ising model, the exponent 6 is 
taken from Refs. |T^. The exponent z in literature ranges from 2.155 to 2.172 
'average' value is given. 
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3 Deterministic dynamics 



Up to now, we have discussed stochastic dynamics induced by Monte Carlo algo- 
rithms. However, in general stochastic dynamics is not equivalent to real dynam- 
ics described by microscopic deterministic equations of motion. Only in some 
cases stochastic dynamics is an effective description of real physical systems. 
Therefore, it is important and interesting to investigate whether there exists uni- 
versal scaling behavior in deterministic dynamics. 

We consider the 2D 0^ theory. The Hamiltonian on a square lattice is 
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and it leads to the equations of motion 

0. = E( 
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(8) 



(9) 



In dynamic evolution governed by Eq. (^, energy is conserved. The solutions 
are assumed to generate a microcanonical ensemble. The temperature could be 
defined as the averaged kinetic energy. In the short-time dynamic approach, 
however, the total energy is a more convenient controlling parameter, since it is 
conserved and can be input from the initial state. Therefore, from now r will be 
understood as a reduced energy density (e — ec)/ec. Here ec is the critical energy 
density. 

The order parameter of the 0^ theory is the magnetization. For the dynamic 
process starting from a random state |]TB|, we assume a same dynamic scaling 
form as in stochastic dynamics. For sufficiently large L and small enough mo and 
t, from scaling form (|l]) one deduces 

(10) 



M{t) = morF(t^/^V 



Similar as the discussion on Eq. (|^) in Sec. 2, we determine the critical energy 
ec and the exponent 9 from Eq. (|l^) and the exponent 1/vz from its derivative. 
Finally, the exponent j3/vz and the dynamic exponent z are estimated from the 
second moment and the auto-correlation. 

In Fig. 1^, numerical solutions of the magnetization with uiq = 0.015 for pa- 
rameters m^ = 2. and A = 0.6 have been plotted with solid lines for three energy 
densities e = 20.7, 21.1 and 21.5 in log-log scale. Careful analysis of the data 
between t = 50 and 500 leads to the critical energy density ec = 21.11(3) [0. 
This agrees well with ec = 21.1 given by the Binder cumulant in equilibrium in 



Ref. |jT6[. At ec, one measures the exponent 9 = 0.146(3). The magnetization 
with mo = 0.009 is also displayed in Fig. ^ (dashed line). The corresponding 
exponent is 6' = 0.158(2). Extrapolating the results to rriQ = 0, we obtain the 
final value 9 = 0.176(7). With the critical energy in hand, we estimate other 
exponents. In Table ^, we summarize all the critical exponents of the 0^ theory. 
Remarkably, not only the static exponents, but also the dynamic exponents z of 
the 0^ theory is the same as those of the Ising model with standard Monte Carlo 
dynamics. 

9 d/z-9 {d-2[3/v)/z l/vz z Ifjjv v 

~lj/ 0.176(7) 0.755(5) 0.819(12) 0.492(26) 2.148(20) 0.24(3) 0.95(5) 
Ising 0.191(1) 0.737(1) 0.817(7) 2.155(3) 1/4 1 

Table 3: The critical exponents of the 0^ theory in comparison with those of the 
Ising model Ref. |H]. 



4 Concluding remarks 

We have demonstrated that universal scaling behavior emerges already in the 
macroscopic short-time regime of the critical dynamic evolution. The short-time 
dynamic scaling is found not only in stochastic dynamics but also in microscopic 
deterministic dynamics. Therefore, it is fundamental. Furthermore, it leads to 
new methods for numerical measurements of both dynamic and static critical 
exponents. The methods do not suffer from critical slowing down. 
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Figure 2: (a)The chiral magnetizations of the 2D FFXY model starting from an 
ordered state. From above, the sohd hnes represent Mi{t) at T=0.452, 0.454 and 
0.456. The dotted hne is at T^ = 0.4545. (b) The derivative (9^/nM/(t, r)|^=o- 
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